Excitonic condensation in the topological insulator bilayer systems: A self-consistent 

mean-field theory study 
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We propose a generic topological insulator bilayer (TIB) system to study the excitonic conden- 
sation with self-consistent mean-field (SCMF) theory. We show that the TIB system presents 
the crossover behavior from the Bardeen-Cooper-Schrieffer (BCS) in weak coupling limit to Bose- 
Einstein condensation (BEC) in the strong coupling limit in zero temperature limit. Moreover, 
comparing to those in the traditional semiconductor systems, we find that the superfluid property 
in the BEC phase is more sensitive to electron-hole density imbalance and the BCS phase is more 
robust in the weak coupling limit. Applying this TIB model into Bi2Se3 material, we find that 
only the BEC phase can be observed in experiment. We also calculate the critical temperature 
for Bi2Se3TIB system and find that is about fOO K. More interestingly, we can expect the relative 
high-temperature excitonic condensation since the Kosterlitz-Thouless transition does not suppress 
the critical temperature for Bi2Se3 in SCMF approximation. 
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I. INTRODUCTION 

Recent technological advances in microfabrication 
bring growing interests in studying exciton condensation 
in different bilayer physical systems such as the semicon- 
ductor electron-hole bilayer [H-Q and graphene bilayer 
0, number of novel physical phenomena are ob- 

tained in these systems, such as the BCS-BEC crossover 
as well as the subtle phase transition in the crossover 
region induced by the density imbalance 0, the dark 
and bright excitonic condensation under spin-orbit cou- 
pling [8|, anomalous exciton condensation in high Lan- 
dau levels in magnetic field Q, room-temperature su- 
perfluidity in graphene bilayer Q, etc. The conven- 
tional electron-hole bilayer is fabricated with semicon- 
ductor such as GaAs/AlGaAs/GaAs. The character of 
the semiconductor electron-hole bilayer systems is that 
the electron and hole bands are quadratic ones with dif- 
ferent effective masses, which means missing particle- 
hole symmetry in these kinds of systems and weak stiff 
phase order. Hence, in semiconductor electron-hole bilay- 
ers, the excitonic condensation needs very low tempera- 
ture. Another better candidate for electron-hole bilayers 
is graphene, which has two-dimensional massless linear 
Dirac-band structure in low energy limit. However, the 
coupling between different Dirac-cone structures in the 
same Brilliouin zone brings flaw to graphene to fabricate 
electron- hole bilayer [loj . 

On the other hand, another growing interest in con- 
densed matter physics is the very recent theoretical pre- 
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FIG. 1: (Color online) Left panel: Schematic structure of 
double-well topological insulators in x-y plane. The external 
gates can independently tune the electron and hole densities. 
Right panel: The linear energy dispersion around the Dirac 
point of the electrons and holes. 



diction [TT|, [T2| and experimental verification of the 
topological insulators (TIs) with strong spin-orbit inter- 
action. Several three-dimensional (3D) solids, such as 
Bii_a;Sba; alloys, Bi2Se3-family crystals, have been iden- 
tified to be strong TIs possessing anomalous band 
structures. The energy scale for the surface states of 
these 3D TIs is dominated by the /c-linear spin-orbit in- 
teraction instead of the fc-quadratic ones even as that in 
semiconductor. Especially, the strong TIs surface has sin- 
gle Dirac-cone band structure which is also different from 
graphene. As a result, it is expected that the excitonic 
condensate of these topological surface states probably 
have new characters. 

Inspired by this expectation, as well as by the re- 
cent experimental access to high-quality TI quantum well 
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films on insulating substrate [19|, in this paper we pro- 
pose a TIB model, a gated double TI layers separated by 
an insulating spacer. Using this TIB model, we theoret- 
ically study the excitonic condensation of these surface 
states. We find that the system also presents BCS-BEC 
crossover along with change of the carriers densities in 
zero temperature limit. However, there are two charac- 
ters different from those of conventional excitonic con- 
densation in semiconductor bilayer systems. The first 
is that the BCS phase of TIB is more robust than that 
of the semiconductor bilayer in weak coupling limit; the 
second is that the superfiuidity of the TIB is more sensi- 
tive to the carriers densities imbalance than that of the 
semiconductor bilayer in strong coupling limit. These 
two characters physically root in the fc-linear band dis- 
persion of the TIB. Moreover, by putting this TIB model 
in Bi2Se3 material, we investigate the excitonic conden- 
sation in the Bi2Se3 bilayer and only find the BEC phase 
occurring due to the values of the parameters of the ma- 
terial. The critical temperature of excitonic condensation 
in Bi2Se3 TIB (~ 100 K) is also calculated in the self- 
consistent mean-field (SCMF) approximation, which is 
found to be higher than that in the traditional semicon- 
ductor electron-hole bilayers. More interestingly, we can 
expect the relative high-temperature excitonic condensa- 
tion since the KT transition does not suppress the critical 
temperature for Bi2Se3 in SCMF approximation. 

II. THE TIB MODEL 

The TIB model is schematically illustrated in the left 
panel of Figure [T] Two TI films are separated by an 
insulating spacer of thickness d, and the electron (hole) 
density can be independently tuned by the external gate 
voltage Vi (V2). The linear dispersions of the TIs around 
Dirac points are cartoonishly depicted in the right panel 
of Figure [TJ A similar setup has been proposed by Zhu et 
al. Q for the conventional semiconductor system. The 
only difference is that the two semiconductor quantum 
wells are now replaced by the two TI films. However, 
this replacement changes the energy dispersion of the 
electrons (holes) and raises different physical properties. 

The grand-canonical Hamiltonian describing this TIB 
system can be written as 

H ^Yl ^'^^'^ ~ 4<TCk^ + ^ (1) 

k,cr 

^ X! ^k-k'C]L+q/2cr'^-k+q/2cr''^-k'+q/2cr'Ck'+q/2(T- 
k,k',q 
a. a' 

Here, k, k', and q are 2D wave vectors in the 
layers, is the quantization volume, cj^^ (cko-) are 
the creation (annihilation) operator for an electron in 
electron (e) layer or hole (h) layer distinguished by 
a-=(e, h). The surface states of the strong TI film 
have the linear dispersion: eke,/i=±fi^WF|k| and /Xg- are 



chemical potentials for electron layer or hole layer. 

Fourier transform of the Coulomb inter- 
action: the intralayer Coulomb repulsive interaction 
V^''{V^'')=2TTe^/ (qe), and the interlayer Coulomb at- 
tractive interaction V^''=— 27re^ exp (— qd) / (qe). Here, 
e is the background dielectric constant and d is width of 
the spacer. 

With the SCMF theory, Eq. ((TJ can be rewritten in a 
2x2 matrix in the basis (e,/i)^, the relevant mean-field 
equations to be solved for the variables /ig, fJ-h, and the 
gap function Ak are 

^'^ = -^E^k-k'^[/(i?k+')-/(i?kO], (2) 

k' ^ 

= ^ E ^k-k' WlfiE+) + vlfiE^,)] , (3) 

k' 

= ^ E ^k-k' HfiK') + ^IfiK')] . (4) 

k' 

= ^ E HfiEi) + vl [fiK)] } , (5) 

k 

"'^ - ^ E {^k[l - fiK)} + vl[l- f{E+)] } , (6) 

k 

where / [E^)^l/{1 + e^^/'=s^) is the Fermi-Dirac dis- 
tribution function, £^k=v^^k+^k '^^^^ ^k=5 {(.ke-^kh), 
Ckcr=ek<T-Mcr+Sk (a'=e, h), and i?^=(5^k±£^k with 
S^k=^ (Cke+Ck/i). In addition, ul=l-vl=^ (l+Ck/£^k)- 
In two dimensions the average interparticle spacing is 
given by [7| 



III. NUMERICAL RESULTS AND 
APPLICATION TO THE BI2SE3 MATERIAL 

Many meaningful physical quantities can be obtained 
by self-consistent solving Eqs. (HJ-®. The exciton's 
energy spectrum is shown in Fig. Ufa) with parameters 
rg=b and a=0 in zero temperature limit. The corre- 
sponding density of states is shown in Fig. HKb). A sta- 
ble energy gap protecting the excitonic condensation is 
evident. Moreover, the wave-vector dependence of Ak 
for equal densities and several values of is shown in 
Fig. 131 We can find the generic feature of the BCS-BEC 
crossover behavior similar to that in the semiconductor 
bilayers. However, the striking character in the TI bilay- 
ers is that the maximum of Ak in the weak coupling limit 
is much stronger than in the traditional semiconductor 
electron- hole bilayers Q. The difference means the BCS 
phase of TIB is more robust than that of the semicon- 
ductor bilayer in weak coupling limit for equal densities 
case. 
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FIG. 2: (Color online) (a) Energy spectrum and (b) density of 
states of the exciton's system with rs=5 and a—Q. We choose 
hvF as the energy unit and nm as the length unit. The diploe 
constant and the spacer width are set as e=d=l. 




FIG. 3: (Color online) Wave-vector dependence of the gap 
function Ak for a=0, d=l, and several values of rg. The 
solid, dashed, dotted, dash-dotted, dash-dot-dotted lines are 
corresponding to rs=0.5, 1.5, 3.0, 5.0, and 20.0, respectively. 
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FIG. 4: (Color online) Maximum value Aniax=max { Ak} as 
a function of a for d=l and several values of ra. The solid, 
dashed, dotted, dash-dotted, dash-dot-dotted lines are corre- 
sponding to rs=0.5, 1.5, 3.0, 5.0 and 20.0, respectively. 



The effect of the density imbalance 
a= {n^-Uh) / [ue+nh) on A^ax is shown in Fig. 
m It is evident to find that density imbalance actually 
suppresses A„iax and that it has different effects on two 
sides of the crossover. In the BEC regime, the main 
effect of the density imbalance reduces the number of 
electron- hole pairs, which results in that the superfluid 
properties are less sensitive to density imbalance. In 
the BCS regime, the density imbalance leads to the 
mismatch of the Fermi surfaces of electrons and holes 
and the finite momentum pairing, which is easier to be 
broken. However, comparing to that in the traditional 
semiconductor bilayers, we find that the superfluid 
property in the BEC phase in our case is more sensitive 
to electron-hole density imbalance. As an example, 
for rs=20 the maximum of gap function Amax for 
TIB disappears at a taking a value smaller than 0.5, 
while it always takes finite values at a varies in the 
whole the zone (—1,1) for the traditional semiconductor 
electron- hole bilayers 0] . 

Now we apply this TIB model to study the conden- 
sation of electron-hole pairs for the topological surface 
states of the Bi2Se3 material. The two TI films in the 
left panel of Fig. [1] now are two ultrathin TI Bi2Se3 
films [13 (about 80 A thick). With the adopted exper- 



imental lattice constants 211 a=4.143 A and c=28.636 



A, we calculate the first-principles surface band struc- 
ture of Bi2Se3 [22] by a simple supercell approach with 
spin-orbit coupling included and obtain the approximate 
Hamiltonian form describing the gapless surface states of 
Bi2Se3 as following, 



(8) 



Although this Hamiltonian has the same form as that of 
the conventional two-dimensional electron gas (2DEG) 
system with Rashba spin-orbit coupling, the intrinsic dif- 
ference between these two kinds of systems is that the 
fc-linear spin-orbit interaction is primary to the TI sur- 
face states, while the parabolic term is dominant in the 
conventional 2DEG. By fitting the first-principles results, 
the parameters in Eq. ([SJ are given as 7=0.21 eV nm^ 
and fivF=0-2 eV nm (namely, w_f=3.04 x lO^m/s). That 
means the energy dispersion around the Dirac point can 
be accurately described by ek=ifiwi?|k| when the wave- 
vector |k| is much smaller than 1.0 nm^^. For numerical 
calculation, we choose nm as the length unit and 0.2 eV 
as the energy unit throughout this paper. The dielectric 
constant e=l and the spacer width d=10 A. In fact, the 
condition that the wave- vector |k| is much smaller than 
1.0 nm^^ requires that only for > 5, then the TIB 
model is valid for Bi2Se3 material. This means that only 
the BEC phase can most probably emerge in the Bi2Se3 
bilayer. 

Now, let us discuss the critical temperature of this 
TIB system. The relation between the maximum value 
Aniax=max {Ak} and temperature T is respectively 
shown in Fig [5] (a) for d=l, a=0, and several values of 
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FIG. 5: (Color online) (a) Maximum value Amax=max{Ak} 
vs temperature T for d—1, a=0, and several values of rs- The 
solid, dashed, dotted, dash-dotted, dash-dot-dotted lines are 
corresponding to rs=0.6, 1.5, 3.0, 5.0, and 20.0, respectively, 
(b) Aniax as a function of the temperature T at rs—5. The 
solid, dashed, dotted, and dash-dotted lines are corresponding 
to a—0, 0.05, 0.1, and 0.2, respectively. 



Ts, and (b) d=l, rs—5.0, and several values of a. From 
Fig. [SJa), we can find that the critical temperature Tc 
decreases as increasing (i.e., decreasing the particle 
density). For the Bi2Se3 bilayer at rs=5.0, the critical 
temperature Tc is calculated as 0.05 in units of 0.2 eV. 
That means the critical temperature Tc is about 8^10 
meV (namely, about 100 K), which is much higher than 
that in the traditional semiconductor electron-hole bi- 
layers. Although the Bi2Se3 TIB system is in the BEC 
phases {rs—5.0, 20.0), the numerical calculated results 
shown in Fig. [SJa) are consistent with the general rela- 
tion of BCS superconductor. 



2A(0) 

Tc 



= 2-Ke 



3.53, 



(9) 



where A(0) is the energy gap at zero temperature. The 
introduced electron-hole density imbalance (a^^O) can re- 
duce the critical temperature. This character is clearly 
shown in Fig. [5Ub): by increasing the density imbalance 
a, the critical temperature Tc decreases. 

As it is known that in two dimensional superfluids, 
the critical temperature is often substantially overesti- 
mated by mean-filed theory. It is ultimately limited by 
entropically driven vortex and antivortex proliferation 
at the Kosterlitz-Thouless (KT) transition temperature 
T'KT=f /5s(2kt) with Ps{T) being the superfluid density 
(the phase stiffness). We can calculate the superfluid 
density by the method of evaluating the couterflow cur- 
rent, which is similar to that in Ref. Joj- The superfluid 
density dependence of temperature is shown in Fig. [6] 
at rs—5 and a=0. From Fig. [51 it is evident to esti- 
mate that the KT transition temperature Tkt is about 
0.05 in units of 0.2 eV. Comparing with the critical tem- 
perature Tc in Fig. [S] at rs=5 and a = 0, the striking 
conclusion is reached: Tc~ Tkt, which means that the 




FIG. 6: (Color online) The calculated Tkt at rs=5 and a=0. 



lOOK high-temperature excitonic condensation may oc- 
cur in the Bi2Se3 TIB system. On the other hand, we can 
estimate the KT temperature with the zero-temperature 
phase stiffness ps{T=0)^Ep /An which is similar with the 
graphene bilayers Q. Considering the case shown in Fig. 
[21 the Fermi energy Ep can be numerical calculated and 
is given as about 0.4 (in units of 0.2 eV). Hence, the KT 
temperature is estimated as TKT~T'i^/8~0.05 in units of 
0.2 eV. Thus, the two estimated methods are consistent 
and the 100 K high-temperature excitonic condensation 
may be protected in the Bi2Se3 TIB system. 



IV. CONCLUSION 



In summary, we perform a generic TIB model to study 
the excitonic condensation with the SCMF theory for 
the topological surface states. Similar to the traditional 
semiconductor electron- hole bilayers, the TIB system 
presents the crossover behavior from the BCS in weak 
coupling limit to Bose-Einstein BEC in the strong cou- 
pling limit in zero temperature limit by changing the ex- 
citon's density. However, two prominent novel characters 
different from the traditional semiconductor electron- 
hole bilayers are found. One is that the superfluid prop- 
erty in the BEC phase is more sensitive to electron-hole 
density imbalance. The other is that the BCS phase 
is more robust than that of the semiconductor bilayer 
in weak coupling limit. Applying this TIB model into 
Bi2Se3 material, we find that only the BEC phase can be 
observed in experiment. Moveover, we theoretically esti- 
mate the critical temperature for the Bi2Se3 TIB system 
and find that it is much higher than that in the tradi- 
tional semiconductor electron-hole bilayers. For example, 
at rs=5 and a=0, the critical temperature Tc is obtained 
as about 100 K. We also study the phase stiffness and 
find that the KT transition does not suppress the critical 
temperature for Bi2Se3 in SCMF approximation. 
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